pub fn new_birthday_probability(n: u32) -> f64 {
    if n >= 365 {
        1.0
    } else {
        let mut res = 1.0;
        let mut i = 1;
        while i < n {
            res *= (365.0 - i as f64) / 365.0;
            i += 1;
        }

        1.0 - res
    }
}
